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SUMMARY 

An integrodif ferential equation that defines the behavior of the shear 
stress in a thermally loaded, thin, flat plate of finite dimensions is derived. 
The equation is reduced to a system of linear algebraic equations by a method 
that employs polynomial approximations for the shear stresses. The boundary 
conditions are satisfied identically. Several examples of the method, presented 
in detail, show that solutions of high accuracy can be produced on a desk calcu- 
lator with a minimum of labor. A calculation of the displacements shows that 
the strain compatibility conditions are closely satisfied everywhere. 


INTRODUCTION 

Several methods now exist for determining the elastic stress distribution 
in a thermally loaded, thin, flat plate of finite dimensions (refs. 1 to 6). 

All these methods are numerical procedures that eventually produce a system of 
simultaneous ordinary differential equations or algebraic equations that require 
considerable labor to effect a solution. 

Reference 6 provides intermediate information from a collocation procedure 
in tabular form that minimizes the effort required to determine the spanwise and 
chordwise stresses for a large variety of plate geometries and temperature dis- 
tributions. However, for span to chord ratios other than those tabulated in 
that reference, an interpolation procedure must be used that requires extreme 
care because the functions vary quite rapidly. A relatively simple method for 
obtaining a solution to this problem is therefore desirable. 

A similar problem .tor a. finite cylinder, rather than for a plate, has 
recently been treated In reference 7. In that reference, an integrodif ferential 
equation in the shear stresses is derived that lends Itself to solution by a 
two-dimensional collocation procedure requiring relatively little labor. 

In a. similar fashion, this report derives the basic integrodif ferential 
equation that defines the behavior of the elastic shear stresses In a finite, 

thin, flat plate subjected to an arbitrary thermal gradient and gives the 
functions necessary to approximate the solution in a collocation procedure. 



Examples that are given and compared with other solutions show good agreement. 
As an additional check, it is shown that the compatibility conditions are 
closely satisfied throughout the region. 


SYMBOLS 

(l 1 ? dimensionless constant 

kZ 

E Young 1 s Modulus 

P^ polynomial in x associated with i^* 1 station 

Q. polynomial in y associated with station 

d 

T temperature function 

T Q reference temperature 

tij constant, equal to dimensionless shear stress at the point (x^yj) 

u spanwise displacement 

v chordwise displacement 

x transformed spanwise coordinate 

y transformed chordwise coordinate 

a linear coefficient of thermal expansion 

J3 span to chord ratio 

j shear strain 

€ normal strain 

7] chordwise coordinate 

v Poisson’s ratio 

| spanwise coordinate 

o normal stress 

T shear stress 



product for all values of k except k = Z 


2 


Superscripts: 

* dimensionless quantity 
' derivative 


ANALYSIS 

Integrodifferential Equation 

With reference to the coordinate system shown in figure l(a)> the relations 
that define the elastic state of stress and strain for the case of plane stress 
are given as follows in reference 8. (Note that the stress-strain relations 
have been modified to include thermal strain terms.) 


= |(°| " V0 T]) + aT 

= ^ a ri - va i ) + aT 


r = 


1 + V 

E 


(lc) 


5T 


a l 



0 


(2a) 




= 0 


ai. 




a 2 


dl 


€_ 
2 ^ 


a 2 


r = o 


Oc 

b 


(2b) 


Making all quantities dimensionless by dividing equations (l) and (3) by 
<xT 0 and equation (2) by EaT Q and transforming I and tj to x and y 

(fig. 1(b)) yield: 


4=4- 

vaj + T* 

(4a) 

e y = °y 

va* + T* 

(4b) 

r* = (i 

+ V )t* 

(4c) 
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1 5 * C) V 

p Sx 0X + T 


= 0 


a * i a 

3y °y + p Sx T * 


= 0 


( 5a) 
(5b) 


a^ e * + a 2 

ay 2 p 2 ax 2 


2 a 2 * 

P ax dy 


= 0 


( 6 ) 


Expressing the compatibility equation, (6), in terms of the stresses, first by 
substituting the stress-strain equations, (4), into equation (6) and then elimi- 
nating the shear terms by differentiating equation (5a) with respect to x and 
(5b) with respect to y and substituting yield: 


l. af_ ag_\ 

p 2 ax 2 ay 2 ) 




(j_ ai_ 

yp 2 a x 2 



( 7 ) 


Integrating equations (5) produces the following expressions for a* and a* 
in the first quadrant: x ^ 


* 

0 x 




dx+ a*(l,y) 


* 

0 


y 




dy+ ctJ(x,1) 


(8a) 


(8b) 


But the boundary conditions on this system require that all normal and shear 
stresses on the free boundaries vanish. Hence, 


Therefore, 


a*(l,y) = a*(x,l) = 0 

* o f S * 

] - t dx 

»x 


( 9 ) 


(10a) 


/ = i f 1 o 

y P 7y £ 


dy 


(10b) 


bubsticuting equations (10) into equation (7) produces a single integrodif i*eren 
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tial equation that defines the shear stresses in the first quadrant of i rj'ure 

1(h): 





T 


* 


( 11 ) 


with the boundary conditions 

T*(x,l) = T*(l,y) = 0 (12) 

Reference G shows that the solution in one quadrant is all that is necessary to 
define the solution everywhere. This is discussed further in the following sec- 
tions . 


Solution by Collocation 

Equation (ll) will be solved approximately by the method of collocation in 
two dimensions. It is assumed that there exists a sequence of functions 
that converge to the solution T* of equation (ll) as m and n approach in- 
finity and that are given by 


m n 

T m (x ' y) = A X) P l U)Q J <y)t ij (13) 

i“l j=l 

where the Pj_ and Qj are known polynomials selected in such a manner that the 
boundary conditions and certain other conditions to be defined are satisfied. 

The constants tqj are as yet unknown. The collocation method now requires 
that the error in replacing t* by in equation (ll) vanish at m by n 

specified points as shown in figure l(b) . To do this, equation (13 ) is substi- 
tuted into equation (ll) and the t— determined so that equation (ll) is sat- 
isfied at each of these m by n stations. This will result in m by n linear 
algebraic equations for the unknown 

In addition to satisfying the boundary conditions the polynomials Pj_ and 
Qj are chosen to satisfy the following conditions: 

P 1 (x i ) = 1 (14a) 

P.(x k ) = 0, k i i (14b) 


The quantity Q,(y) is similarly a polynomial in y associated with the U* 1 
J 

station and satisfies similar conditions: 
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(15a) 

(15b) 


V y j> " 1 

«j(jk) - °> k t 0 

Hence^ as seen from equation (13), the boundary conditions, equations (12), are 
identically satisfied, and tjj = T*(xj_,yj)* 

Polynomials that have the desired properties are easily obtained. For 
example. 



where TT is the product for all values of k except k = i. Polynomials of 
kf i 

the form of equation (16) are of the lowest possible degree that will satisfy 
conditions ( 15 ) . 

Equation (13) is now substituted into equation (ll) and evaluated at each of 
the m stations in x and n stations in y to produce a system of m by n 
simultaneous linear algebraic equations in tj_j of the following form: 

m n 

y t ■ • 

L~i kZ ij 

i=l j=l 

where 

c ki - - 1 p i«j + 3 p i"/ 1 «j d y + wj"/ p i dx < 18 > 

P Jy •'X 

y - n 

The shear stresses at any point can now be determined from equation (13). The 
spanwise and chordwise stresses are given by 

m n 

c{ = ( e E «: 

i=l j=l 3 
m n ( rl 

°y = §■ £ £ P i / Qj <3y tij (19b) 

P i=l j=l Jy 

obtained by substitution of equation (13) into equations (10) . 




?&• 


k = 1 , 2 , . . ,, jn (17) 


2 = 1 , 2 , . 
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Special Cases of Temperature Variation 


If the temperature function is an even or odd function of only one of the 
independent variables, the solution need be determined in only one quadrant of 
the coordinate system shown in figure l(b). The solution everywhere else is 
known by symmetry, or asymmetry , depending on the function. Writing the equa- 
tions in only one quadrant reduces the number of stations required for a given 
degree of accuracy, significantly reducing the amount of labor required to pro- 
duce a solution. Inasmuch as any function of a single variable can be separated 
into the sum of an even function and an odd function, it is possible to solve 
the problem of an arbitrary temperature distribution in one variable in two 
steps and add the two solutions. For the case in which the temperature function 
depends on two variables, it is frequently possible to separate the function 
into the sum of two functions of a single variable. This procedure is covered 
more fully in reference 6. Separating the problem into several simpler problems 
will always result in a reduction in labor. For example, if it is decided to 
use four stations along each axis to represent the solution, 16 simultaneous 
equations of the form of equation (17) will be produced, and 16 terms will ap- 
pear in each double summation. But there are only two stations along each axis 
in each quadrant that produce only four simultaneous equations and only four 
terms in each double sum. Hence, even if as many as four solutions were re- 
quired, an even and odd solution in each variable, less labor is required. An 
example of an even polynomial that has the desired properties is 



An example of an odd polynomial that has the desired properties is 



EXAMPLES 

Two by Two Collocation on a Square Plate 

As a first example, consider a square plate subjected to a thermal gradient 
in the y-direction independent of x: 


T* = y 2 

This temperature distribution will produce a normal stress distribution symmet- 
rical about both axes and a shear stress distribution that is asymmetrical about 
both axes. Therefore, polynomials of the form of equation (21) will be chosen. 
Two equally spaced stations in each direction will be used; and, therefore, the 
P polynomials are identical to the Q polynomials. Writing the polynomials at 
equally spaced points, in this case at l/4 and 3/4, produces 
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P x = Q x = 8.533 r 5 - 13.33 r 3 + 4.800 r 

P 2 = Q 2 = -6.095 r 5 + 6.476 r 5 - 0.3810 r 

where r is a dummy variable representing x for the P. and y for the 
polynomials . Substituting these polynomials into equations (17) and (18) pro- 
duces the following system: 


—45.92 tpp 
96.70 t u 
96.70 tpp 
-60.05 tpp 


- 17.99 t 12 

- 72.17 t 12 
+ 92.10 t 12 
+ 56 • 72 t-^2 


- 17 .99 t 2 p 
+ 92 .10 t 2 p 

— 72.17 t 2 p 
+ 56 • 72 t^p 


+ 12 .37 t 2 g 

- 68.11 t 22 

- 68 • 11 tg 2 

- 64.50 t 22 


- 2.0 


— Ci * 


o 


- 2.0 


- 2.0 


The solution to the system is : 


t 1± = 0.02991 
t 12 = 0.04683 
t 21 = 0.04683 


t 22 


0 .08553 


Substituting these values and the P’s and Q's evaluated at the desired x’s 
and y's into equations (13) and (19) produces the stresses everywhere. The 
stress distribution in the first quadrant is tabulated in columns 4, 7 and 10 
of table I and is shown in figure 2. 

Compared with the distribution produced by equations (13) and (19) of this 
report is the elastic stress distribution in unpublished NASA data given in 
columns 3, 6, and 9 of table I and shown in figure 2. These stresses are cal- 
culated from Airy’s stress function, which is produced as the solution to the 
biharmonic equation. The solution was obtained by a finite-difference method, 
which required the solution of 441 simultaneous linear algebraic equations 0 An 
inspection of table I and figure 2 discloses that the agreement between the two 
methods is excellent. 
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Three by Two Collocation on a Square Plate 


As a second example y consider the square plate of the first example sub- 
jected to the same one -dimensional thermal gradient. For this example^ however , 
take three stations in the x-direction. The P polynomials become 


P ± = -41.66 x 7 + 81.00 x 5 - 46.58 x 3 + 7.232 x 
P 2 = 27.00 x 7 - 46.50 x 5 + 20.02 x 3 - 0.5208 x 
P 3 = 15.90 x 7 + 20.32 x 5 - 4.529 x 3 + 0.1104 x 


written at the equally spaced stations x^ = l/6, x^ = l/2, = 5/6 . Substi- 

tuting into equations (17) and (18), solving the resultant 6 by 6 matrix, and 
evaluating the stresses by using equations (13) and (19) produce the solution 
tabulated in columns 5, 8, and 11 of table I and shown in figure 2. It is 
evident that the additional stations produced no improvement in the solution. 


Two by Two and Three by Two Collocation on a Rectangular Plate 

Next, consider a plate that has a span three times its chord and is subject- 
ed to the same thermal gradient as the square plate of the preceding discussion. 
Using the 2 by 2 station collocation produces the solution tabulated in columns 
4, 7, and 10 of table II and shown in figure 3. Comparing that solution with the 
finite-difference solution tabulated in columns 3, 6, and 9 of table II and shown 
in figure 3, shows poor agreement. However, if another station is chosen in the 
x-direction, the solution is considerably improved. This is also shown in fig- 
ure 3 and columns 5, 8, and 11 of table II. 


Comparison of Displacements 

As an additional check, the compatibility of strains was determined from the 
strain-displacement relations. Relations were derived for both spanwise and 
chordwise displacements in two different ways. First, the appropriate strain- 
displacement equation was integrated directly (e.g., eq. (22b) for the spanwise 
displacement). Then the shear -displacement equation (22c) was integrated with 
fespect to the appropriate variable, and the normal displacement was eliminated 
by making use of the remaining equation (e.g., for the spanwise displacment eq. 
(22c) was integrated over y, and the chordwise displacement was eliminated by 
using eq. (22a)). The displacements were then calculated by the two different 
relations and compared. 

Proceeding in this manner, from reference 8, 
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(22a) 


* 1 du 

x p ax 

1 / du 1 dv* 

2 Idy + P dx 


Integrating equation (22b) yields 


u? = 3 


* , 
€_ dx 


The spanvise displacement may be calculated a different way by using equations 
(22a) and (22c): 


* o 

u 2 = 2 


'0 JQ 


r 

/ 

Jo 


dy + f^x) 


From equation (22a), 


* / * , 

v = / e y dy 


Therefore, 


/.y /.y 

i tl 

3 / J 

•/n •'o 


u! = 2 / T* dy 


/ 0 •'o 


dy dy + f^x) 


or integrating by parts yields 


<-* /V*-J f £ 


dy + 


if 


y &* 

y dy + f n (x) 


vhere 


fn (X) = U*(X,0) 


To make ut = % at y = 0, f (x) is chosen to equal u*(x,0). Hence, 
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( 23 ) 


* f J * 

u n = 3 / e x dx 

*'0 


* 

= 2 


/ y /*y v * 

r " dy -f/ lr 

Jo 


sr dy + f 


At* 

Jo 


dy + uq(x,0) 


(24) 


Similarly, 


v i = 


/ 


y * 

€ y dy 


(25) 


v* = 2p / T* dx- 3 2 x / dx + 3 2 

''O "V '0 


/ 


del 


f 


de* * 

x dx + v 1 (0^y) 


(26) 


For a parabolic thermal gradient, T* = y c and v = l/3, from equations (4), 
(13), and (19) the result is 


e* J 1 

y 


' m n f 1 

E Ep Qj 

f i=i j=i L J 


o m n 

- 1 E E Qj 

1=1 j=i 




: 1 tij + Y c 


’ 27a) 


e x = 



m n 


33 


E E P: 

i=l j=l 



Qj dy + y c 


(27b ; 


F = | EnEVu 


(27c; 


Substituting equations (27) into equations (23), (24), (25), and (26) produces a 
check on the compatibility of the strains. Figure 4 and columns 3 to 6 of 
table III show values of u*, u*, v*, and v* throughout the first quadrant for 
2 by 2 collocation on the square plate . It is evident that there is excellent 
agreement between the displacements. 


DISCUSSION 

The examples of the preceding section demonstrate the simplicity and accu- 
racy of this method. Reference 6 points out that for a one -dimensional tempera- 
ture function, it is unnecessary to use values of 3 greater than 3, which 
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produces the solution for the semi -infinite plate at the center. For most prob- 
lems, therefore, it is unnecessary to use more than three stations in the x- 
direction, indicating that the largest system of equations to be solved will be 
6 by 6. This permits the use of a desk calculator instead of an electronic data 
processing machine, as required by most of the other methods discussed* 

Despite the fact that the entire stress distribution is written in terms of 
the shear stresses at only four or six points within the region of definition 
and those stresses are calculated from an equation that is approximated only at 
those same points, a solution of remarkable accuracy is obtained* At this 
point, it is appropriate to mention that previous uses of the double collocation 
method (ref* 9) required more stations, because the solution was obtained in 
Airy's stress function, which required the solution of a higher order equation 
and also required differentiation to obtain the stresses. The solutions pre- 
sented in references 5 and 6 are compared with those in references 1 to 4 and 
with each other. Therefore, it was necessary to compare the solutions produced 
herein with only one of the others. The finite-difference solution was chosen 
because the only approximations used there were the standard ones used in numer- 
ical quadrature and differentiation. It is to be noted that the solutions prod- 
uced by the various references and the solution of this report are virtually 
identical. Where any differences occur, the check of the compatibility condi- 
tion produced by the method reported herein places greatest confidence in this 
solution . 

Because of the small number of stations used in the approximation, the tem- 
perature function must not change rapidly. This restriction was also pointed 
out in reference 6. However, it is impossible to state quantitatively how rap- 
idly the function may vary. For the case of a rapidly varying function, more 
stations must be taken. One way to assure that a sufficient number of stations 
has been chosen is to solve the problem several times for increasing numbers of 
stations. When increasing the number of stations produces no change in the so- 
lutions, it can be assumed that maximum accuracy has been obtained from the 
method. References 5 and 6 state that for the problems considered, up to a 
ninth-power temperature function, no improvement was noted for more than three 
stations in one direction. It should be noted that in those references, col- 
location was in one direction only, producing a system of simultaneous linear 
ordinary differential equations in Airy's stress function. 


CONCLUSIONS 

A method of determining the elastic state of stress in a thermally loaded^ 
thin, flat plate of finite dimensions has been presented, which, for the class 
of temperature distributions investigated, produces a very accurate solution 
with a minimum of labor. Its principal advantage is that it produces a solution 
of the same order of accuracy of other methods by a technique that requires the 
solution of only four or six simultaneous linear algebraic equations instead of 
a system of many algebraic equations or a system of differential equations. Its 
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limitations are that the temperature function may not be a rapidly varying 
function and must be separable into the sum of functions of a single variable. 


Levis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, April 8, 1963 
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TABLE I. 


DIMENSIONLESS STRESS IN SQUARE PLATE SUBJECTED TO CHORDWISE PARABOLIC THERMAL GRADIENT 




Coorc 

iinate 



Stress distribution 



Shear stress. 


Span, 

Chord, 

Snanwise . 


nVi -? <3 0 



T* 


X 

y 


4 



* ' 


Finite- 

rJ1 f* f* 

2 by 2 

3 by 2 



Finite - 
difference 
solution 

2 by 2 
colloca- 
tion 

3 by 2 
colloca- 
tion 

Finite - 
difference 
solution 

2 by 2 
colloca- 
tion 

3 by 2 
colloca- 
tion 

solution 

colloca- 

tion 

colloca- 

tion 






Column 






1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

0 

0 

0.141 

0.141 

0.141 

0.141 

0.141 

0.142 

0 

0 

0 


.2 

.131 

.131 

.131 

.131 

.130 

.131 

0 

0 

o 


.4 

.095 

.096 

.096 

.102 

.101 

.102 

0 

0 

o 


.6 

.016 

.017 

.017 

.061 

o 060 

.060 

0 

0 

0 


.8 

-.136 

-.138 

-.138 

.020 

.020 

.020 

0 

0 

0 


1.0 

-.409 

-.412 

-.411 

0 

0 

0 

0 

0 

0 

0.2 

0 

0.131 

0.130 

0.130 

0.131 

0.131 

0.131 

0 

0 

o 


.2 

.122 

.122 

.121 

.122 

.122 

o 122 

.020 

.020 

.020 


.4 

.089 

.090 

.090 

.095 

.095 

.095 

.035 

.036 

.036 


.6 

© 01 6 

.017 

.017 

.057 

.057 

.057 

.042 

.043 

.043 


• 8 

-.126 

-.128 

-.128 

.019 

.019 

.019 

.034 

.034 

.034 


1.0 

-.384 

-.386 

-.386 

0 

0 

0 

0 

0 

0 

0.4 

0 

0.102 

0.101 

0.101 

0.095 

0.096 

0.095 

0 

0 

o 


.2 

.095 

.095 

.095 

.089 

.090 

.089 

.035 

.036 

.036 


.4 

.071 

.072 

.072 

.071 

.072 

.072 

.064 

.065 

.065 


.6 

• 015 

.016 

.016 

.044 

.045’ 

.045 

.078 

.079 

.079 


o 8 

-.099 

-.100 

-.100 

.016 

.016 

.016 

.063 

i064 

.064 


1.0 

-.311 

-.312 

-.312 

0 

0 

0 

0 

0 

0 

CD 

O 

0 

0.061 

0.060 

0.060 

0.016 

0.017 

0.016 

0 

0 

o 


.2 

.057 

.057 

.057 

.016 

.017 

.016 

.042 

.043 

.042 


.4 

.044 

.045 

.045 

.015 

.016 

.015 

.078 

.079 

.078 


.6 

.011 

.012 

.012 

.011 

.012 

.011 

.096 

.098 

.098 


.8 

-.059 

-.060 

-.060 

.005 

.005 

.005 

.081 

.082 

.082 


1.0 

-.196 

-.197 

-.196 

0 

0 

0 

0 

0 

0 

0.8 

0 

0.020 

0.020 

0.020 

-0.136 

-0.138 

-0.137 

0 

0 

o 


.2 

.019 

.019 

.019 

-.126 

-.128 

-.127 

.034 

.034 

.034 


.4 

.016 

.016 

.015 

-.099 

-.100 

-.100 

.063 

.064 

.064 


.6 

• 005 

.005 

.005 

-.060 

-.060 

-.061 

.081 

.082 

.082 


.8 

-.021 

-.020 

-.020 

-.021 

-.020 

-.020 

.069 

.071 

.071 


1.0 

- • 065 

-.069 

-.069 

0 

0 

0 

0 

0 

0 

1.0 

0 

0 

0 

0 

-0.409 

-0.412 

-0.412 

0 

0 

0 


.2 

0 

0 

0 

' -.384 

-.386 

-.387 

0 

0 

0 


.4 

0 

0 

0 

-.311 

-•312 

-.311 

0 

0 

0 


.6 

0 

0 

0 

- © 196 

-.197 

-.195 

0 

0 

0 


.8 

0 

0 

0 

- © 065 

-.069 

-.068 

0 

0 

0 

1 

1.0 

0 

0 1 

0 

0 

0 

0 

0 

0 

0 



TABLE II 


DIMENSIONLESS STRESS IN RECTANGULAR PLATE SUBJECTED TO CHORDWISE PARABOLIC THERMAL GRADIENT 


Coordinate 

Stress distribution 

Shear stress, 

T* 

Span, 

X 

Chord, 

y 

Spanwise, 

a x 

Chordwise, 

Finite- 

difference 

solution 

2 by 2 
colloca- 
tion 

3 by 2 
colloca- 
tion 
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Dimensionless chord coordinate, y 
(b) Chordwise stress. 

ire 2. - Concluded, Stress in square plate subjected to 
Lordvise parabolic thermal gradient. 
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Dimensionless span coordinate, x 
(a.) Spanwise displacement. 

Figure 4. - Displacements in square plate subjected to chord 
wise parabolic thermal gradient. 
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Dimensionless chord coordinate, y 
(b) Chordwise displacement. 

Figure 4. - Concluded. Displacements in square plate subjected to chord 
wise parabolic thermal gradient. 
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